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The dark energy component of the universe is often interpreted either in terms of a cosmological 
constant or as a scalar field. A generic feature of the scalar field models is that the equation of 
state parameter w = P/p for the dark energy need not satisfy w = — 1 and, in general, it can be 
a function of time. Using the Markov chain Monte Carlo method we perform a critical analysis of 
the cosmological parameter space, allowing for a varying w. We use constraints on w(z) from the 
observations of high redshift supernovae (SN), the WMAP observations of CMB anisotropics and 
abundance of rich clusters of galaxies. For models with a constant w, the ACDM model is allowed 
with a probability of about 6% by the SN observations while it is allowed with a probability of 
98.9% by WMAP observations. The ACDM model is allowed even within the context of models 
with variable w: WMAP observations allow it with a probability of 99.1% whereas SN data allows it 
with 23% probability. The SN data, on its own, favors phantom like equation of state (w < —1) and 
high values for Qnr- It does not distinguish between constant w (with w < —1) models and those 
with varying w(z) in a statistically significant manner. The SN data allows a very wide range for 
variation of dark energy density, e.g., a variation by factor ten in the dark energy density between 
z — and z = 1 is allowed at 95% confidence level. WMAP observations provide a better constraint 
and the corresponding allowed variation is less than a factor of three. Allowing for variation in w has 
an impact on the values for other cosmological parameters in that the allowed range often becomes 
larger. There is significant tension between SN and WMAP observations; the best fit model for one 
is often ruled out by the other at a very high confidence limit. Hence results based on only one 
of these can lead to unreliable conclusions. Given the divergence in models favored by individual 
observations, and the fact that the best fit models are ruled out in the combined analysis, there is 
a distinct possibility of the existence of systematic errors which are not understood. 
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I. INTRODUCTION 



Observational evidence for accelerated expansion in 
the universe has been growing in the last two decades 
. Observations of high redshift supernovae 0, 01 pro- 
vided an independent confirmation. Using these along 
with observations of cosmic microwave background ra- 
diation (CMB) 0, and large scale structure [(J Q, 
we can construct a "concordance" model for cosmology 
and study variations around it (e.g., see [a, JH ; for an 
overview of our current understanding, see |lO|). 

Observations indicate that dark energy should have an 
equation of state parameter w = P/p < —1/3 for the 
universe to undergo accelerated expansion. Indeed, ob- 
servations show that dark energy is the dominant com- 
ponent of our universe. The cosmological constant is the 
simplest explanation for accelerated expansion 0, 0] 
and it is known to be consistent with observations. In 
order to avoid theoretical problems related to cosmologi- 
cal constant i other scenarios have been investigated. 
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In these models one can have w — 1 and in general w 
varies with redshift. These models include quintessence 
[l3| . k-essence 0], tachyons ^1^3, phantom fields [l7| . 
branes ITa . etc. There are also some phenomenological 
models |l9f . field theoretical and renormalisation group 
based models (see e. g. J20l| 'l . models that unify dark mat- 
ter and dark energy |2l[ and many others like those based 
on horizon thermodynamics (e.g. see |22|). Even though 
these models have been proposed to overcome the fine 
tuning problem for cosmological constant, most of these 
require similar fine tuning of parameter(s) to be consis- 
tent with observations. Nevertheless, they raise the pos- 
sibility ofw(z) evolving with time (or being different from 
— 1), which — in principle — can be tested by observa- 
tions. 

Given that w < —1/3 for dark energy for the universe 
to undergo accelerated expansion, the energy density of 
this component changes at a much slower rate than that 
of matter and radiation. Indeed, w = — 1 for cosmological 
constant and in this case the energy density is a constant. 
Unless w is a rapidly varying function of redshift and be- 
comes w ~ at (z < 1), the energy density of the dark 
energy component should be negligible at high redshifts 
(z 3> 1) as compared to that of non-relativistic matter. If 
dark energy evolves in a manner such that its energy den- 
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sity is comparable to, or greater than the matter density 
in the universe at high rcdshifts then the basic structure 
of the cosmological model needs to be modified. We do 
not consider such models here. We confine our attention 
to models with dark energy density being an insignificant 
component of the universe at z > 1 and choose observa- 
tions which are sensitive to evolution of w(z) at redshifts 

To put the present work in context, we recall that com- 
bining supernova observations with the WMAP data pro- 
vides strong constraints on the variation of dark energy 
density [23] • (A review of relevant observations for con- 
straining dark energy models along with a summary of 
the previous work in this area is given in Section Iff Bl ) 
Reproducing the location of acoustic features requires the 
angular diameter distance to the last scattering surface to 
be in the correct range. This analysis showed that while 
the data from SN observations allows for a large range in 
parameters of dark energy, combining with WMAP data 
limits this range significantly. However, in that work, we 
did not explore the cosmological parameter space widely 
and had fixed nearly all parameters other than those used 
to describe evolution of dark energy. In the present work, 
we allow many cosmological parameters to vary and in- 
clude constraints from cluster abundance in addition to 
the supernova and WMAP constraints. 

In addition to obtaining quantitative bounds on pa- 
rameters in different contexts, we address the following 
key issues in this paper: 

• Does allowing cosmological parameters to vary 
weaken the constraints on variation of dark energy? 

• Conversely, how does the allowed range for different 
cosmological parameters change when we allow for 
a epoch dependent w(z) ? 

• Do the observational constraints agree with each 
other? In particular, what kind of cosmological 
models are preferred by SN and WMAP observa- 
tions individually ? 

The last point is important and requires elaboration. 
Different observational sets are combined together pre- 
cisely because these observations are sensitive to differ- 
ent combinations of cosmological parameters and facili- 
tate in breaking degeneracies between parameters. If we 
consider ACDM models then the SN observations, for ex- 
ample, broadly depend on the combination (0.85^1 nr~ 
0.531V) while WMAP is sensitive to (Sl NR +n v ) ||, 
a feature which was originally highlighted in the litera- 
ture as 'cosmic complementarity'. Therefore, we cannot 
expect constraints from different observations to agree 
over the entire parameter space. At the same time, we 
do not expect models favored by one observation to be 
ruled out by another when such a divergence is not ex- 
pected. This divergence may point to some shortcomings 
in the model, or to systematic errors in observations, or 
even to an incorrect choice of priors. If all observational 



sets are consistent then we should be able to derive simi- 
lar constraints using subsets of observations, even though 
the final constraints may not be as tight as with the full 
set of observations. 

In order to address the questions listed above in a sys- 
tematic manner, we proceed in three steps. We choose a 
'base' reference model with cold dark matter and cosmo- 
logical constant, with neutrinos contributing a negligible 
amount to the energy density of the Universe. We as- 
sume that the Universe is flat and restrict ourselves to an 
unbroken power law for the primordial power spectrum 
of density fluctuations and we assume that the pertur- 
bations are adiabatic. Another assumption is that the 
perturbations in tensor mode are negligible and we take 
r = [2^]. We choose this to be our standard model as 
this can be described by a compact set of parameters. 

Next, we generalize from ACDM models (w = —I) to 
study a wider class of dark energy models with a con- 
stant w and address the issues listed above. In this case, 
we also study the effect of perturbations in dark energy. 
Finally, we generalize to models in which w is allowed to 
vary with z in a parameterized form. This approach al- 
lows us to delineate changes that come about from choos- 
ing a constant w =/= —I, from those allowed by a varying 
dark energy. We do not impose theoretically motivated 
constraints on models, e.g. we do not require w > — 1 as 
the present work is focused on understanding the nature 
of models favored by observations. 

The paper is organized as follows: In section[n]we dis- 
cuss the background cosmological equations followed by a 
brief review on the various observation used to constrain 
dark energy equation of state and the observations we 
concentrate on. The Markov Chain Monte Carlo method 
is discussed in section IIIII and detailed results are pre- 
sented in section llVl We conclude with a discussion of 
the results and future prospects for constraining dark en- 
ergy models in section [V] 

II. THEORETICAL BACKGROUND 
A. Cosmological equations 

If we assume that each of the constituents of the ho- 
mogeneous and isotropic universe can be considered to 
be an ideal fluid, and that the space is flat, the Friedman 
equations become: 

iS) - 

d AttG , „ 

- = ^(P + 3P) (2) 

a 3 

where P is the pressure and p = Pnr + p 7 + p DE with 
the respective terms denoting energy densities for nonrel- 
ativistic matter, for radiation/relativistic matter and for 
dark energy. Pressure is zero for the non-relativistic com- 
ponent, whereas radiation and relativistic matter have 
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P 7 = p 1 /'d. If the cosmological constant is the source of 
acceleration then p DE = constant and P DE = — p DE . 

An obvious generalization is to consider models with 
a constant equation of state parameter w = P j p = con- 
stant. One can, in fact, further generalize to models with 
a varying equation of state parameter w{z). Since a func- 
tion is equivalent to an infinite set of numbers (defined 
e.g. by a Taylor-Laurent series coefficients), it is clearly 
not possible to constrain the form of an arbitrary func- 
tion w(z) using finite number of observations. One pos- 
sible way of circumventing this issue is to parameterize 
the function w(z) by a finite number of parameters and 
try to constrain these parameters by observations. There 
have been many attempts to describe var y in g dark ene rgy 
with different parameterizations pi Ea El El El l3c| 
where the functional form of w(z) is fixed and the vari- 
ation is described with a small number of parameters. 
Observational constraints depend on the specific param- 
eterization chosen, but it should be possible to glean some 
parameterization independent results from the analysis. 

To model varying dark energy we use two parameteri- 
zations 



)(z) = u> + w'(z = 0) 



(1 + z)p 



p=l, 2 (3) 



These are chosen so that, among other things, the high 
redshift behavior is completely different in these two pa- 
rameterizations [23] . If p = 1 [3l| , the asymptotic value 
w(oo) = wq + w'(z = 0) and for p = 2, w(oo) = Wq. 
For both p = 1,2, the present value w(0) = wq. Clearly 
we must have w(z S> 1) < —1/3 for the standard cos- 
mological models with a hot big bang to be valid. This 
restriction is imposed over and above the priors used in 
our study. 

Integrating d(pa 3 ) = —w(z)pda , the energy density 
can be expressed as 



Pde _ /j _|_ z yi(l+w +w' a ) 



Pde 

for p = 1 (in Eq. and 



Pde 



(1+*) 



3(l+i«o) 



exp 



exp 



3u>n 



3w' Q )z 
' 1 + z 



(4) 



1 + z 



(5) 



for p — 2. The allowed range of parameters Wq and 
Wq = w'(z = 0) is likely to be different for different p. 
However, the allowed variation at low redshifts in p DE 
should be similar in both models as observations actu- 
ally probe the variation of dark energy density. Indeed, 
in an earlier study |23j where we had studied a restricted 
class of models, we found this to be the case. For exam- 
ple, p DE can vary by at most a factor two up to z = 2 
when both the WMAP and SN data are taken into ac- 
count p3| . This reaffirms the expectation that the results 
are parameterization independent at some level. 



B. Observational constraints 

In this subsection, we briefly review potential observa- 
tional constraints on dark energy and we also summarize 
previous work in this area. 

Constraints on dark energy models essentially arise 
as follows: To begin with, dark energy affects the rate 
of expansion of the universe and thus the luminosity 
distance and also the angular diameter distance. Con- 
straints from observations of high redshift supernovae 
H |H HI |H H HI and the location of peaks in the 
angular power spectrum of CMB anisotropics mainly use 
this feature j^ij • The signature of acoustic peaks in cor- 
relation function of galaxies also provides a similar geo- 
metric constraint (37]]. There is also an effect on weak 
lcnsing statistics through changes in distance-redshift re- 
lation]!! m m m Ea . 

Second, the rate of expansion influences the growth of 
perturbations in the universe and this leads to another 
set of probes of dark energy ^3|- Abundance of rich 
clusters of galaxies, their evolution and the integrated 
Sachs- Wolfe (ISW) effect belong to this category of con- 
straints, along with constraints from redshift space dis- 
tortions UlllillilliJIilli^. theS<3 constraints are 
sensitive to different aspects of dark energy and a com- 
bination of all of these can put tight limits on models. 
The redshift space distortions are a local effect as these 
are sensitive to the rate of growth of density perturba- 
tions at a given epoch. The abundance of rich clusters 
of galaxies, the ISW effect and distances are integrated 
effects in that the effect of dark energy is averaged over 
a range of redshifts in some sense. 

Observations of high redshift supernovae of type la 
provide the most unambi guo us evidence for accelerated 
expansion 0, 0, 133, El l34j. Assuming these sources 
to be standard candles, observations spanning a range of 
redshifts can be used to study the change in rate of expan- 
sion and this imposes direct constraints on the variation 
of dark energy density. Supernovae have been observed 
up to a redshift of z max ~ 1.8 and hence can be used 
to constrain models of dark energy up to this redshift. 
Constraints from SN observations alone however, permit 
a large variation in the dark energy parameters and 
in particular favor models with w < — 1 at the present 
epoch 0,0 EH- 

Baryon oscillations in the matter-radiation fluid prior 
to decoupling provide a standard scale and the angle at 
which the acoustic peaks occur in the angular power spec- 
trum of temperature anisotropies in the CMB fixes the 
distance to the surface of last scattering. This provides a 
useful constraint on models of dark energy |3(j as long as 
dark energy does not affect the dynamics of universe at 
the time of decoupling of matter and radiation. Unlike 
supernovae that are observed over a range of redshifts 
where dark energy is dominant, the surface of last scat- 
tering is at z ~ 1100. However, the exquisite quality of 
CMB anisotropy measurements makes this a very useful 
constraint and these observations offer a constraint that 
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is different from SN observations [2j| . Indeed, as we shall 
see, WMAP and SN observations often favor models that 
are mutually unacceptable. 

Recent detection of the baryon acoustic peak in galaxy 
correlation function using the luminous red galaxy sam- 
ple of the SDSS survey has provided an additional handle 
to constrain cosmological parameters [37j- The geomet- 
ric constraint from these observations can, in principle, 
constrain models of dark energy. A measurement of the 
angular scale corresponding to the peak at different red- 
shifts can indeed be a powerful constraint. 

If we consider a given cosmological model that is con- 
sistent with observations of CMB anisotropy then the 
amplitude of fluctuations at the time of decoupling is 
fixed, and its linearly extrapolated value today can be 
computed using linear perturbation theory. The abun- 
dance of rich clusters of galaxies is related to the ampli- 
tude of perturbations in dark matter at a scale of about 
8h _1 Mpc. If we study different models for dark energy 
while other parameters are not changed, the abundance 
of rich clusters constrains the net growth of structures 
between the epoch of decoupling and the present epoch 

Redshift space distortions due to kinematics and the 
Alco ck- Paczy nski effect are also potential probes of dark 
energy IfH] • Ongoing surveys like the SDSS and fu- 
ture surveys will be able to distinguish between differ- 
ent dark energy models through these effects How- 
ever, this method does not provide useful constraints at 
present. 

Dependence of the distance-redshift relation and a dif- 
ferent rate of growth for perturbations, as well as changes 
in the matter power spectrum are also reflected in weak 
lensing statistics. Several studies have been carried out 
on the potential constraints that can be put on dark en- 
ergy models from weak lensing observations and their de- 
generacy with other parameters. These studies indicate 
that future surveys will be able to put strong constraints 
on dark energy models [3^. El) . 

Growth of perturbations also leaves a signature in the 
CMB anisotropy spectrum at large angular scales. The 
ISW effect leads to an enhancement in the angular power 
spectrum at these scales. The detailed form of this en- 
hancement depends on the equation of state parameter 
w and its variation. This effect can be detected by cross- 
correlation of temperature anisotropies with the fore- 
ground distribution of matter |53j . It is difficult to distin- 
guish the ISW effect from the effect of a small but non- 
zero optical depth r due to re-ionisation by using only the 
temperature anisotropies in the CMB, cross-correlation 
with the matter distribution or polarization anisotropics 
in the CMB must be used. 

Redshift surveys of galaxy clustering do not constrain 
properties of dark energy directly, however the shape 
of the power spectrum constrains the combination T = 
^NRh |54|. This provides an indirect constraint on dark 
energy through the well known degeneracy between £Inr 
and w (e.g. see p^l. 



The large number of different observations that can 
be used to constrain dark energy models is encouraging. 
Indeed, many attempts have been made to use some of 
these observations to put constraints on models 0, I23I 



C. A choice of three observations 

In this work, we concentrate on SN, WMAP and clus- 
ter abundance observations. We briefly explain the rea- 
son for this choice and the kind of constraints one can 
expect. 

The left panel in Fig. 2] shows the degeneracy in fl nr 
and wo (for models with constant w; i.e., with w = 
in Eq. |3J). The figure shows contours of constant lumi- 
nosity distance Hodi(z) at z = 0.17 (red/solid curves) 
and z = 1.17 (blue/dashed curves). The second panel 
displays the constant luminosity distance contours in the 
w'(z — 0) — wq plane for p = 2 if £Inr = 0.3. Given 
that SN observations constrain luminosity distance as a 
function of redshifts, these figures illustrate the shape 
of the allowed region that we are likely to get and also 
demonstrates the degeneracies between different parame- 
ters. The third panel shows how the redshift at which the 
expansion of the universe begins to accelerate depends on 
the parameters wo and w'(z = 0) for p = 2. This epoch 
is constrained by SN observations and hence the allowed 
region in parameter space should lie between contours of 
this nature. Clearly, regions with a late onset of acceler- 
ation (upper right corner) as well as a very early onset 
of acceleration (lower right corner) will be ruled out by 
observations of supernovae. 

The fourth panel of this figure shows the redshift at 
which matter and dark energy contribute equally in terms 
of the energy density of the universe. Structure formation 
constraints are likely to follow these contours as the rate 
of growth for density perturbations is significant only in 
the matter dominated era. Too little structure formation 
(upper right corner) as well as too much structure forma- 
tion (lower left corner) are likely to restrict the allowed 
models along a diagonal (upper left to lower right) in this 
plane. 

Lastly, the location of acoustic peaks in the angular 
power spectrum of temperature anisotropies in the CMB 
is the most significant constraint provided by CMB ob- 
servations. This essentially constrains the distance to 
the surface of last scattering and hence a suitably de- 
fined (see eqn.JSJ) mean value (w e ff) for w. The right 
panel shows contours of w e f / , which run almost diagonal 
in this plane. Thus a band of allowed models is the likely 
outcome of comparison with observations. The contours 
of w e ff are the same as contours of equal distance to 
the surface of last scattering, or the I corresponding to 
the first peak in the angular power spectrum of CMB 
temperature fluctuations. 

SN data provides geometric constraints for dark energy 
evolution. These constraints are obtained by comparing 



FIG. 1: The left figure in this panel shows contours of constant luminosity distance HoDl(z)/c in £Inr — w plane. The 
red/solid contours are for luminosity distance at redshift z = 0.17 and blue/dashed ones are for luminosity distance at redshift 
z = 1.17. The red/solid contours from top right to bottom left correspond to values between HqDl{z = 0.17)/c = 0.19 — 0.21 
and blue/dashed lines are for HoDl(z = 1.17)/c = 1.52 — 2.45. The other plots show contours in the wo —w(z = 0) plane. The 
second figure shows contours of constant HoDl(z = 0.17)/c and HoDl(z = 1.17)/c. Starting from upper right the contours 
are for HoDl(z — 0.17)/c = 0.18 — 0.2 and HoDl(z = 1.17)/c = 1.63 — 2.1. The third panel shows contours of redshift at 
which the expansion of the universe starts to accelerate. Starting from upper right corner in counter-clockwise direction, the 
contours are for z = 0, 0.2, 0.4, 0.6 and 0.8. The constraint from SN observations essentially constrains this redshift and hence 
the allowed region can be expected to follow these contours. The next panel shows red-shift at which the energy density in 
matter and dark energy are equal. The curves, starting from left are for z — 0.2, 0.4, 0.6, 0.8 and 1.0. Structure formation 
constraints are likely to follow these contours as the rate of growth for density perturbations is significant only in the matter 
dominated era. The right panel shows contours of w e ff at the surface of last scattering (see Eq. From top right towards 
bottom left the values are w e ff = —0.5, —1.0, —1.5 and —2.0 If the location of acoustic peaks provides the main constraint 
then the region allowed by CMB observations should follow these contours. We have used fijvfl = 0.3 for these plots. 



the predicted luminosity distance to the SN with the ob- 
served one. The theoretical model and observations are 
compared for luminosity measured in magnitudes: 



m B (z) = M + 5log 10 (D L ) (6) 



where M = M - 5log 10 (H ) and D L = H d L , M be- 
ing the absolute magnitude of the object and is the 
luminosity distance 



f dz 

d L = (1 + z)a{t )r(z); r(z) = c j (7) 



where z is the redshift. This depends on evolution of 
dark energy through H(z). For our analysis we use the 
combined gold and silver SN data set in 0. This data 
is a collection of supernova observations from 0, H3 and 
many other sources with 16 supernovae discovered with 
Hubble space telescope Q • The parameter space for com- 
parison of models with SN observations is small and we 
do a dense sampling of the parameter space. 

CMB anisotropics constrain dark energy in two ways, 
through the distance to the last scattering surface and 
through the ISW effect. Given that the physics of recom- 
bination and evolution of perturbations does not change 
if w(z) remains within some safe limits, any change in 
the location of peaks will be due to dark energy (3(|. 
The angular size 9 of the Hubble radius at the time of 



decoupling can be written as: 

e _ x = H(z)/H 

}dy/(H(y)/H ) 
o 

\/n NB ,(i + zf 

fdy/Jn NR (l + z) 3 + B DE (z)/Q$ E 



/ dy/yJfl NR (l + zf + n de (1 + z f 1+w '" ] 



The second line is obtained as decoupling happens at a 
redshift where dark energy is not important, and if we 
ignore the contribution of radiation and relativistic mat- 
ter; the last equation defines W e ff. Clearly, the value of 
the integral will be different if we change wq, w'(z = 0) 
and there will also be some dependence on the param- 
eterized form. Location of peaks in the angular power 
spectrum of the CMB provide a constraint, but this can 
only constrain w e ff and not all of wq, w'(z — 0) and 
p. Therefore if the present value wq < Weff then it is 
essential that w'(z = 0) > 0, and similarly if wq > w e tf 
then w'(z — 0) < is needed to ensure that the integrals 
match. Specifically, the combination of wo, w'(z = 0) 
and p should give us a w e ff within the allowed range. 

In our analysis, we use the angular power spectrum of 
the CMB temperature anisotropies [i(J E3 as observed 
by WMAP and these are compared to theoretical pre- 
dictions using the likelihood program provided by the 
WMAP team [5^ ■ We vary the amplitude of the spec- 
trum till we get the best fit with WMAP observations. 
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TABLE I: This table lists the priors used in the present work. 
Apart from the range of parameters listed in the table, we 
assumed that the universe is flat. We assumed that the pri- 
mordial power spectrum had a constant index. Further, we 
ignored the effect of tensor perturbations. The range of val- 
ues for wo and w'iz — 0) is as given below, but with the 
constraint that w(z — 1000) < —1/3. Any combination of 
wo and w'(z — 0) that did not satisfy this constraint was not 
considered. 



Parameter 


Lower limit 


Upper limit 


n B 


0.03 


0.06 




0.1 


0.5 


h 


0.6 


0.8 


r 


0.0 


0.4 


n 


0.86 


1.10 


w 


-2.0 


-0.4 


w'(z = 0) 


-5.0 


5.0 



III. MARKOV CHAIN MONTE CARLO 
METHOD 

We compute y 2 using the routines provided by the 
WMAP team The CMBFAST package [H is used 
for computing the theoretical angular power spectrum for 
a given set of cosmological parameters. We have com- 
bined the likelihood program with the CMBFAST code 
and this required a few minor changes in the CMBFAST 
driver routine. Given the large number of parameters, 
the task of finding the minimum x 2 and mapping its be- 
havior in the entire range of values for parameters is com- 
putationally intensive. 

We adapt the Metropolis algorithm 60] (also known in 
the context of parameter estimation as the Markov Chain 
Monte Carlo (MCMC) method |H E2) for efficiently 
mapping regions with low values of % 2 . The algorithm 
used is as follows: 



Note that this is different from the commonly used ap- 
proach of normalizing the angular power spectrum at 
I = 10. As we use the entire angular power spectrum 
for comparison with observations, the impact of ISW ef- 
fect on the likelihood is relatively small. 

It has been pointed out that constraints from struc- 
ture formation restrict the allowed variation of dark en- 
ergy in a significant manner I5U . We use observed abun- 
dance of rich clusters |44ll45ll46| to apply this constraint. 
Since the mass of a typical rich cluster corresponds to the 
scale of 8h _1 Mpc, cluster abundance observations there- 
fore constrain erg, the rms fluctuations in density contrast 
at 8h _1 Mpc. The number density of clusters depends 
strongly on <jg and Qnr- We use the cr 8 constraints given 
in 44] from ROSAT deep cluster survey and are given 
by cr 8 = (0.58 ±0.1) x n N ° R i7+0 - imNR at 99% confidence 
level. The cosmological model should predict erg in the 
allowed range in order to be consistent with observations. 

Recent detection of the baryon acoustic peak using 
luminous red galaxy sample of the SDSS survey has 
provided an additional handle to constrain the cosmo- 
logical parameters j^j. We also used distance scale 
Dy at redshift z = 0.35 introduced in the above refer- 
ence to further constrain the cosmological models. Here 
Dy(0.35) 3 = Dm(z) 2 cz/H(z) and the observational con- 
straint is £V(0.35) = 1370±64 Mpc at the la level. This 
fourth observation does not add significantly to other 
constraints listed here and we will not describe quan- 
titative results from this constraint here. 

Priors used in the present study are listed in Table 1. 
Apart from these limits on the models studied here, we 
also assumed that neutrinos are massless and the ratio 
r of tensor to scalar mode is zero. These assumptions 
are consistent with the known upper bounds, and in any 
case these do not make any difference to the observations 
used here as constraints I25L 



1. Start from a random point in parameter space 
and compute C; and x 2 ( r i)- 

2. Consider a small random displacement r.; + i = V{ + 
dr and compute x 2 ( r i+i)- 

3. If X 2 ( r i+i) < X 2 ( r i) then i — > i + 1. Go to the first 
step. 

4. Else: 

• Compute Ax 2 = x 2 ( r i+i) — X 2 ( r i) an d 
exp[— a Ax 2 ]- 

• Compare this with a random number < (3 < 
1. 

• If P < exp[— a Ax 2 ] then i — * i + 1. Go to the 
first step. 

The size of the small displacement dr and the parame- 
ter a are chosen to optimally map the regions of low x 2 - 
We wish the chain to converge towards the minimum, 
starting from an arbitrary point, and we also want the 
Markov chain to map the region in parameter with low 
X 2 exhaustively without getting bogged down near the 
minima. These two conflicting requirements are recon- 
ciled by choosing a small but non-zero value of a. Maxi- 
mum displacement allowed in one step is small compared 
with the range of parameters, but not small enough for 
the chain to get trapped in a small region around the 
minimum. The optimum values of maximum displace- 
ment and a are related to each other. We ran several 
chains with a varying number of points, a typical chain 
has about 10 4 points. For each set, we have at least 10 5 
points (We have done calculations for five sets: cosmolog- 
ical constant, constant w with and without perturbations 
in dark energy, time varying w(z) for p — 1 and p = 2. 
Results presented here required an aggregate CPU time 
of nearly 10000 hours on 2.4 GHz Xeon CPUs). The 
convergence criteria for such chains is satisfied for all the 
sets, and for all the parameters in each set |fj2l ]. 
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We use the MCMC approach only for comparison of 
models with the CMB data. Observations of cluster 
abundance are compared with models from the Markov 
chain run for CMB, after the chain has been run. Com- 
parison of models with observations of high redshift su- 
pernovae is done separately. This approach is more con- 
ducive to one of the questions that we wish to address, 
namely, are the observational constraints consistent with 
each other? 



IV. RESULTS 

We present results in the form of likelihood functions 
for various parameters in sets of increasing complexity, 
starting with the standard ACDM model. Before we pro- 
ceed with a discussion of results in this form, we discuss a 
few specific models sampling a few interesting regions of 
the parameter space in order to develop an intuitive feel 
for different observational constraints. We call these fidu- 
cial or reference models. Along with the fiducial models, 
we also discuss the best fit models in each set. We find 
the best fit model for individual observations as well as 
for the combination of all the observations. 



A. Fiducial Models 

1. The ACDM model 

The ACDM model is our 'standard' model and we be- 
gin our discussion with this class of models. Several stud- 
ies have been carried out to constrain parameters for the 
ACDM model HQ. Our results for the ACDM model 
bring out — among other things — the differences from 
previous work which arises due to a different method we 
use here for normalizing power spectra. (See section 2.3 
for details.) Differences introduced by priors are also ap- 
parent. Our results are as follows: 

• For ACDM model, if we consider SN observations 
alone, we get a best fit at FInr = 0.28 with a 
X% = 233.1. (We will use subscript S for x 2 from 
SN analysis and W for analysis with WMAP ob- 
servations.) This model with FInr — 0.28 is al- 
lowed by WMAP observations and has a best fit 
Xw = 974.3 for tt B = 0.045, h = 0.69, n = 0.95 
and r = 0.008. SN observations do not fix these 
parameters so we varied the other parameters to 
get the best fit WMAP model for fl NR ■ = 0.28. 

• The model which best fits the WMAP observations 
has n B = 0.05, Q NR = 0.34, h = 0.66, n = 0.96 
and r = 0.002 with a, xlv = 972 - 5 - The X 2 value 
corresponding to SN fit is Xs ~ 239.9. In the 
context of cosmological constant models alone, this 
model is away from the SN best fit by = 6.8 
and is allowed with probability 0.009. In other 



words, the model most favored by WMAP observa- 
tions is allowed by the SN observations only with 
less than one percent probability (We define prob- 
ability V of a given model to be 1 — C/100, where 
C% is the confidence limit at which the model is 
allowed. By using this definition we avoid dilu- 
tion due to a large parameter space. While the 
statement about x 2 is accurate and directly obtain- 
able from the analysis, the conversion of confidence 
intervals to probabilities has well-known statisti- 
cal caveats while dealing with multiparameter fits. 
This should be kept in mind while interpreting our 
statements about probability with which a model 
is allowed). In contrast, the model most favored by 
SN observations is allowed by WMAP observations 
with a probability V = 0.945. 

We now restrict some of the parameters to values fa- 
vored by other observations e.g. j^. We fix the baryon 
density parameter fl R = 0.05, present day Hubble pa- 
rameter h — 0.7 and the spectral index n = 1 for 
these models. Allowing CInr and r to vary the best fit 
ACDM model in this restricted class of models is with 
fl NR = 0.31 and r = 0.14 and Xw = 974 - 8 - This is 
fairly close to the best fit model found by the WMAP 
team using a large set of observations This model is 
allowed by the rich cluster abundance observations, and 
also by SN observations (xg — 234.8, the corresponding 
probability being V — 0.2). 

Thus convergence between the WMAP and SN obser- 
vations happens only in a narrow window for flat ACDM 
models, with the SN constraint being the tighter of the 
two. It is worth mentioning that in a wider class of mod- 
els, (obtained by relaxing the prior fitot = 1) SN data 
favors a closed universe with il to t — 1-44 ± 0.28 and — 
more importantly — allows the f2 tot = 1 models with 
V = 0.12 

2. Models with a constant w 

We now allow the dark energy equation of state pa- 
rameter to have values different from w = — 1 but do not 
allow variation with time. We then find that: 

• SN observations favor a model with w = —1.99 
and tt NR = 0.47 and with x% = 227.5. This is the 
root cause of the phantom menace. This model is, 
however, ruled out at a very high probability by 
WMAP data (with a Axw = 13.6) and is allowed 
only with V = 0.022. 

• WMAP observations favor higher values for w 
(non-phantom models) and the best fit model has 
w = —0.72. The other cosmological parameters 
corresponding to this best fit are £l R = 0.06, 
n NR = 0.34, h = 0.64, n = 0.99 and r = 0.19 
with Xw = 971.6. SN observations, on the other 
hand, rule out this model at a very high significance 
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TABLE II: This table lists best fit parameter values and \ 2 f° r different models and some of the selected fiducial models. The 
abbreviation b.f. denotes the best fit model for the particular data set. The corresponding probabilities are given in brackets. 
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level (Ax| — 53). This is the case when we do not 
take perturbations in dark energy into account for 
computing theoretical predictions. 

• If we include dark energy perturbations [31 the 
model which fits best with WMAP observations 
is same as the best fit for the ACDM models. 
This model is allowed by SN observations with 
X 2 = 239.3 which is allowed with probability 0.003. 
Clearly, the discrepancy between the WMAP and 
SN observations is reduced when we take pertur- 
bations in dark energy into account, but SN ob- 
servations allow the WMAP best fit model only 
marginally even in this case. 

• How well does the ACDM fare when we allow a 
range of w ? The cosmological constant is allowed 
with V = 0.063 by SN observations in the set of 
models with constant w. WMAP observations al- 
low the ACDM models with a very high probability, 
indeed the best fit model continues to be a ACDM 
model if perturbations in dark energy are taken into 
account. 

To illustrate these aspects, we consider two fiducial 
models with w = —1.5 and w — —0.9 allowing £Inr 
and t to vary. Other parameters are fixed as for the 
ACDM model (see subsection 1 above). We do not take 
perturbations in dark energy into account here. The best 
fit values using WMAP observations for w = —1.5 are 
Qnr — 0.38 and r = 0.04. Generically, models with 
lower w require higher values of Qnr which is clear from 
the form of the curves in Fig^](a). This model is allowed 
by WMAP observations with Xw = 978.8 as well as by 
SN observations (x% = 229.9). But this model is not 
allowed by observations of cluster abundance. The value 
of erg — 1-1 f° r this model is higher than the upper limit 



erg = 1.02 allowed with 99% confidence level for this value 

Of Q, NR . 

With w = -0.9, WMAP favors a model with Q NR = 
0.28 and r = 0.18 {x 2 w = 974.5). This model is allowed 
by cluster abundance observations. The model becomes 
marginal when SN observations are taken into account 
with x% = 237.9, as compared to a minimum of Xs = 
227.5 for models with constant w. 

Once again, we find that WMAP observations and 
SN observations favor different regions in the parameter 
space. Generalizing to the class of models with a constant 
w from w = — 1, wc find that SN data has a distinct pref- 
erence for w < — 1 and the ACDM model is allowed only 
marginally. WMAP data also shows a mild preference 
for models with w ^ — 1, though it continues to allow 
the ACDM model. If perturbations in dark energy are 
taken into account then the ACDM model continues to 
be the most favored model for WMAP observations. 



3. Models with varying w(z) 

Next, we allow the dark energy equation of state pa- 
rameter to vary. As it is not possible to take the effect of 
non-adiabatic perturbations into account for these mod- 
els and as we do not wish to add another parameter, 
we work with models without any perturbations in dark 
energy [64j. It is clear that this will introduce a slight 
bias towards models with w > — 1 but will not change 
anything else. For the discussion of fiducial models, we 
choose the parameterization with p = 2 in Eq. Q . 

• Parameter values for the best fit model with SN 
data are wq = —1.95, w'(z — 0) — —4.5 and 
fl NR = 0.498 with a x| = 227.4. As with the 
constant w models, SN data favors large negative 
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values of the equation of state parameter. SN ob- 
servations do not favor models with varying w(z) 
over models with w ^ — 1, since the change in best 
fit x 2 is only 0.1 when the additional parameters are 
added. WMAP observations allow this model with 

V = 0.035 (xw = 987.6) if we have n B = 0.056, 
h = 0.662, n = 1.04 and r = 0.002. However, this 
model is ruled out by cluster abundance observa- 
tions as eg = 1-25 for this model is higher than the 
allowed range at 99% confidence level. 

• The best fit model for WMAP data has w Q = 
-1.48, w'(z = 0) = 3.86, fl B = 0.05, Q NR = 0.24, 
h = 0.73, n = 1.1 and r = 0.35 with Xw = 970.9. 
This model is allowed by SN observations with 

V = 0.33 (Ax 2 = 4.63). 

The tension between the WMAP and SN observations is 
less serious for models with varying w than for models 
with a constant w. Part of the reason is that with a 
larger number of parameters, the same Ax?; gives us a 
larger probability for a given model. 

As in the previous cases, let us choose fiducial models 
by restricting some of the parameters. We consider a 
model with dark energy parameter values wq = —1.5 
and w'(z = 0) = 1.0. The WMAP best fit with these 
parameters is with £Inr = 0.37 and r = 0.0005 with 
Xw = 976.5 . The change in acceptance level of this 
model is due to our restricting the values of the spectral 
index and the Hubble parameter. This model has Xs = 
229.3 and is allowed by all the three observations used 
here. 

If we move closer to the model most favored by SN 
observations, w n = —1.5 and w'(z = 0) = —5.0, then we 
find that WMAP data favors Q NR = 0.42 and r = 0.0. 
Even this model is outside the range allowed by WMAP 
observations at 68% confidence limit and is allowed at 
V = 0.12. The model fares similarly with SN observa- 
tions, i.e., it is allowed with V = 0.12. 

We wrap up with a discussion of a model with wq > 
— 1. We consider u> = —0.9 and w'(z = 0) = —3.0. 
WMAP observations allow this model with Xw = 977.4 
for Qnb, — 0.365 and r = 0.0013. SN observations al- 
low this model within the 95% confidence limit with a 
probability V — 0.15. This model is also allowed by ob- 
servations of cluster abundance. 

Within the context of models with a variable w(z), 
the ACDM model is allowed by WMAP observations 
(Axw = 1.6) as well as by supernova observations 
(Ax| = 5.73). SN observations clearly favor models 
other than the ACDM model in context of U, tot = 1, while 
no such preference is seen for WMAP observations. How- 
ever, models favored by SN observations require £Inr. to 
be much larger than the values favored by observations 
of rich clusters 0]. 

In summary, there is significant tension between the 
sets of observations we are studying and this tension does 
not reduce when the parameter space is enlarged. We 
see that best fit model for one set of observation is often 



ruled out by another set with a high level of significance. 
There is an overlap of region allowed with 95% confidence 
limit in all cases and within 68% confidence limit in some 
cases; but this does not take away the significance of the 
differences which the above analysis has thrown up. The 
results presented here are summarized in Table 2. 

B. Results in Detail 

In this section we outline the detailed results of our 
analysis. We marginalize the results in the multi- 
dimensional parameter space to derive likelihood func- 
tion for each parameter of interest. The likelihood func- 
tion is sensitive to the bin-size used for the given param- 
eter and tends to be somewhat noisy. It is customary 
to smooth the likelihood function with a Gaussian filter 
in order to remove noise but the results are sensitive to 
the width of the filter used. Therefore we choose to plot 
the cumulative likelihood as this is insensitive to binning 
and smoothing is not required. We define the cumulative 
likelihood as follows: 

T ( x ) = S ( 9 ) 

where x is the parameter we are interested in and it has a 
range x m m < x < x m ax, £(x) is the likelihood obtained 
by marginalizing over other parameters. The central 
value for the given variable is thus x c , with I(x c ) = 0.5. 
This can be, and in general it is different from the value 
of the parameter for maximum likelihood. I(x) is like the 
cumulative probability function for the parameter x. Pa- 
rameter values for which I(x) = 0.025 and I(x) = 0.975 
define the range allowed at 95% confidence limit, i.e. the 
probability that the variable lies within this range is 0.95. 
We mark this limit by two vertical lines in all the likeli- 
hood plots. 

C. The ACDM model 

Fig. |21 shows marginalized cumulative likelihood I(x) 
for the ACDM model, different frames correspond to the 
different parameters we have considered here. Curves are 
shown for Qnr, the shape parameter T — £Inr h, Hub- 
ble parameter h = H /100 kms^Mpc" 1 , fi& h 2 , spec- 
tral index n and optical depth to the redshift of reioni- 
sation r. Red/dot-dashed curves show constraints from 
WMAP observations of CMB temperature anisotropies, 
dark-green/dotted curve shows constraints from WMAP 
observations and abundance of rich clusters, black/solid 
curve shows constraints from a combined analysis of 
WMAP observations, cluster abundance and high red- 
shift supernovae. The blue/dashed curve shows the con- 
straints from SN observations alone, this has been plot- 
ted only for Vl^r as SN data does not constrain other 



10 




FIG. 2: This Figure shows plots of I (see eqnQJJ for various parameters for the ACDM model. Red/dot-dashed curves show 
constraints from the WMAP observations of temperature anisotropies, dark-green/dotted curve shows constraints from WMAP 
observations and abundance of rich clusters, black/solid curve shows constraints from a combined analysis of WMAP obser- 
vations, cluster abundance and high redshift supernovae. The blue/dashed curve shows the constraints from SN observations 
alone, this has been plotted only for the relevant parameters. Vertical lines mark the 95% confidence limit for each parameter. 



parameters directly. Combined analysis for other param- 
eters does have an input from supernova observations 
as many models allowed by WMAP and cluster abun- 
dance observations are ruled out by SN observations. 
Vertical lines mark the 95% confidence limit when all 
the observations are used together. CMB observations 
allow considerable range for each of these parameters 
whereas observations of high redshift supernovae provide 
tight constraints. The reason why SN observations pro- 
vide a tight constraint is clear from the discussion in the 
previous section, namely, SN observations favor mod- 
els with n, tot > 1. Within the context of models with 
^tot = 1, SN observations favor models with w <C — 1. 
The ACDM model is only marginally allowed by SN ob- 
servations within both the sets; flat ACDM models are 
allowed with V = 0.12 |24| and within flat models the 
ACDM model is allowed with V = 0.06. 

The allowed range for all the parameters within 68% 
confidence limit is given in Table 3 and range allowed in 
95% confidence limit is given in Table 4. Table 3 clearly 
illustrates the discrepancy between SN observations and 



WMAP observations. The values are comparable with 
those obtained in other analyses p| The range of 
values for the shape parameter T — CInrIi favored by 
these observations is consistent with values obtained from 
galaxy surveys 0, IU E3 • The allowed range of values for 
h are in agreement with direct determination |63| . 

Constraints from abundance of rich clusters do not 
make a significant impact on the likelihood function of 
individual parameters even though these constraints re- 
ject a significant fraction of models allowed by WMAP 
observations. To illustrate this, we have plotted points 
in the Markov chain for ACDM models that are allowed 
by WMAP observations in a few projections in the pa- 
rameter space (Fig.[3J). Also shown in the same plots are 
points allowed by abundance of rich clusters of galaxies. 
This clearly shows that the region in parameter space 
allowed by cluster abundance is distinctly smaller than 
that allowed by CMB observations. Well known degen- 
eracies between parameters are also highlighted by this 
figure, e.g., there is a clear degeneracy between £Inr and 
h. An important point to note is that if the preferred 
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FIG. 3: This Figure shows points in the Markov chain for ACDM models that are allowed by WMAP observations of temperature 
anisotropies in the CMB within 95% confidence limit in the full parameter space (pink points). Models that are also allowed 
within a 99% confidence limit by abundance of rich clusters of galaxies are shown as dark green points. The points are 
shown in a few projections to illustrate the fact that the requirement of cluster abundance rejects many models allowed by 
CMB observations. While for some projections it is clear that the region in parameter space allowed by cluster abundance is 
distinctly smaller than that allowed by CMB observations, in other projections these seem to have an almost complete overlap. 
This plot also highlights degeneracies in the parameter space. 



range for aa were to be towards larger values than taken 
here 0, |6g then the cluster abundance constraint will 
favor models with larger Qnr- There is also a related 
shift towards w < — 1. 



D. Models with a constant w 

We now consider models with a constant equation of 
state parameter w. Introduction of this additional pa- 
rameter changes the relative effectiveness of different ob- 
servations in constraining cosmological parameters. Ob- 
servations of high redshift supernovae of type la constrain 
the key cosmological parameters much more strongly 
than CMB observations for ACDM models. This is no 
longer the case once we introduce w as a parameter. The 
main reason for this is the degeneracy between w and 
Qnr- 

In models with w ^ 1, it is necessary to take per- 
turbations in the dark energy component into account 
mi l67ll68| . Here, we study models with and without per- 
turbations in dark energy in order to illustrate the role 
played by these perturbations and to study how strongly 
these influence determination of cosmological parame- 
ters. Fig. 0] shows the likelihood / for parameters Qnr, 
w, h, flf, h 2 , n and r in models with and without per- 
turbations in the dark energy component. Models with 
a larger Qnr. and smaller w are better fits to SN obser- 
vations, whereas CMB observations prefer models with 
smaller Qnr and a larger equation of state parameter w. 
The combination of these observations and abundance of 
rich clusters constrains both the parameters to a fairly 
narrow range, much narrower than is allowed by SN ob- 
servations alone. Models with w > — 1 fare badly with 



CMB observations when perturbations in the dark en- 
ergy component are taken into account. For models with 
a constant w, observations allow higher values of Qb h 2 
than for ACDM models. The allowed range in r is smaller 
then the case with no dark energy perturbations whereas 
the ranges for n are similar in both. 

We find that the range of w allowed at the 95% confi- 
dence limit is smaller when perturbations in dark energy 
are allowed. For other parameters, the allowed range is 
similar. In other words, if we ignore perturbations in 
dark energy we can still make a reasonable estimate of 
the range of parameters allowed by observations. This 
fact is of immense use when we work with models that 
have a varying equation of state parameter w. In or- 
der to take full effects of dark energy perturbations in 
these models it is essential to know full details of the 
model 0, |63|. We cannot include the effect of non- 
adiabatic perturbations in a model independent study of 
dark energy models with a varying w. Given the fact 
that ignoring perturbations in dark energy does not lead 
to an incorrect estimate of the range of parameters (ex- 
cept for w) that is allowed, we can safely proceed with 
our analysis without taking perturbations in dark energy 
into account. As regards the equation of state parameter 
w, ignoring perturbations tends to allow w > — 1 models 
with a larger probability and we should keep this in mind 
while interpreting results. 

E. Models with varying w(z) 

We now proceed to the case of varying w, we use two 
parameterizations given in |23| . The first of these, corre- 
sponding to p — 1, is a Taylor series expansion for w in 
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FIG. 4: This figure shows I for Unr, w and h for cosmological models with a constant equation of state parameter w. Panels 
in the first and third row are for models without perturbations in the dark energy component, whereas panels in second and 
fourth row are for models with perturbations in dark energy. Each panel shows / derived from observations of high redshift 
supernovae, WMAP observations, WMAP observations combined with the constraints from abundance of rich clusters, and all 
three observations in combination (with the same color coding as in Fig. 2). We have marked 95% confidence limits derived 
from using all the observations in concert as two vertical lines. The relevant parameters are written in the x-axis labels of the 
figures. 
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FIG. 5: This figure shows / for parameters that describe the 
variation of the equation of state parameter: wo and w'(z = 
0). Left panels are for p = 1 and the right ones are for p — 2. 
Color coding used is same as in figure 2. 
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FIG. 6: This figure shows I for w e ff at the redshift of surface 
of last scattering for the models with constant w (blue/dashed 
curve), variable w with p — 1 (green/dot-dashed curve) and 
with p = 2 (purple/dotted curve). Perturbations in dark 
energy are not taken into account in these models. See text 
for a more detailed description of this figure. 



scale factor and this is a very commonly used parameter- 
ization. The variation of the equation of state parameter 
is monotonic in this case and rapid increase of w at low 
redshifts cannot be allowed as it will lead to w > —1/3 at 
high redshifts. The parameterization with p = 2 avoids 
this problem to some extent as the value of w at very 
high redshifts is the same as the present value, but there 
can be a large deviation from this at low redshifts with 
the deviation peaking at z = 1. Rapid variation at low 
red-shifts has been reported [34| on the basis of SN ob- 
servations 0, , and even though these conclusions have 
been contested 69] it is useful to check if the larger set 
of observations support a rapid variation of w at low red- 
shifts. 

We have already seen in the discussion of fiducial mod- 
els, supernova observations do not distinguish between 
models with a constant w and models with a variable 
equation of state parameter. Also, WMAP observations 
do not differentiate between the three classes of models 
being studied here in a statistically significant manner. 

Fig. shows the allowed ranges of parameters wq and 
w'{z = 0) for p = 1 (left panels) and p — 2 (right panels). 
In both the parameterizations, the preferred values for 
u>o with supernova observations are for phantom models 
(wq < — 1) and a tendency for a larger w at intermediate 
redshift, i.e., w' > though supernova observations do 
not provide a clear constraint on this parameter. This is 
perhaps related to the fact that there is a strong degener- 
acy in Qnr and wq. CMB observations and abundance of 
rich clusters of galaxies allow models around the ACDM 
model, which is fairly close to the centre of the allowed 
region. A combination of these three observations re- 
jects models with wq <C — 1 due to CMB constraints and 
wo > — 1 due to SN constraints. Even though all the 
observations allow w' = 0, the combination of these ob- 
servations does not favor such models. This implies that 
the overlap of allowed regions for the three observations 
is stronger for models with w' ^ 0, it is clear from table 3 
that there is no overlap between SN and WMAP at 68% 
confidence limit for constant w. The ACDM model, i.e., 
wq = —1 and w = is a marginally allowed model for 
both the parameterizations. 

To understand the nature of constraint from CMB ob- 
servations [hJ], we computed the likelihood of w e ff (as 
defined in eqn.JSJ) for models by comparing these with 
WMAP data. This is then compared with the likelihood 
for w in models with a constant equation of state pa- 
rameter. We have plotted this in Fig. [S] for p = 1 as 
well as p = 2, along with the curve for constant w (with 
no perturbations in dark energy). All the three curves 
show very similar behavior and the 95% confidence limit 
is identical for all three (—1.5 < w e ff < —0.6). This 
also shows that the CMB observations primarily provide 
a constraint for w e ff. Given that adding perturbations 
reduces the likelihood for models with w > —1, it is likely 
that detailed analysis of a model with perturbations in 
dark energy taken into account will limit the range for 
w e ff in this region. 
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FIG. 7: Panels for £Inr, h, T = Q.nr h, Qt h 2 , n and r. Red/solid curve is for ACDM, blue/dashed curve for models with 
constant w, green/dot-dashed for p = 1 and purple/dotted curve for p = 2. 



Lastly, we study the effect of varying dark energy on 
other parameters. The specific question we wish to ad- 
dress is, how the allowed ranges for these parameters 
change if we allow variation of dark energy. Fig. shows 
likelihood for the parameters studied here. We have plot- 
ted the likelihood using all three observations for the 
ACDM models, constant w models as well as for p = 1 
andp = 2. For most parameters, the effect of w ^ —1 and 
varying dark energy is to increase the range of allowed 
values. This increase in the allowed range is sometimes 
accompanied by a shift, e.g. for h where varying dark 
energy models fit observations better with smaller values 
as compared to the ACDM model as well as models with 
constant w. This shift is primarily due to models with 
w > — 1 and this point has been noted in other analyses 
as well . If this is the case then including perturbations 
in dark energy may well remove this shift. 

Similarly, larger values of spectral index n and opti- 
cal depth to the epoch of reionisation r fit observations 
better. As these parameters can be constrained using 
other observations, we may be able to restrict models 
with varying w by constraining the values of these pa- 
rameters. For example, polarization anisotropies in the 
CMB can be used to constrain t mEll. We find that 



the presently available information from WMAP about 
polarization anisotropies does not lead to a significantly 
improved constraint on the parameter r. 



F. Evolution of dark energy 

We now summarize our results for the allowed vari- 
ation of dark energy, once all three observational con- 
straints have been taken into account. In the left panel 
of Fig. we have plotted the cumulative likelihood for 
the equation of state parameter w at redshift z = 1. 
Here we have used models which lie within the range 
—1.1 < wq < —0.9. The upper and lower panels corre- 
spond to parameterizations with p = 1 and p = 2 respec- 
tively. The allowed range of variation in the equation 
of state by supernova observations is much larger than 
that allowed by WMAP results. This, again, is a reflec- 
tion of the strong preference of supernova observations 
for w -C — 1 and of the large parameter space allowed 
by SN data. (Our result that SN data prefers w <C — 1 
with large variation is consistent with previous published 
analysis e.g., in H3-) The range of values in both the pa- 
rameterizations are similar for this subset of models. In 
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FIG. 8: The top and bottom panels in the figure correspond to parameterizations with p = 1 and p = 2 respectively. The 
figures on the left show likelihood for the equation of state parameter w at redshift 2 = 1 for models which have equation of 
state lying within the range —1.1 < wo < —0.9. In the middle panel we plot the likelihood for variation in the equation of state 
parameter from the present value at redshift z = 1. The third panel shows the allowed range of change in the dark energy 
density upto = 1. Color coding is same as in Fig. 2 



the middle panel we have shown the likelihood for varia- 
tion in the equation of state parameter from the present 
to its value at redshift z = 1 . The allowed ranges of vari- 
ation in dark energy equations of state are different for 
these two parameterizations. In fact, the constant dark 
energy equation of state is ruled out at 95% confidence 
level for p = 2 (that is, the probability of occurrence is 
less than 0.05) when all the constraints are taken into ac- 
count even though each observational constraint allows 
such models individually. Clearly, the models with con- 
stant w allowed by each of these observations are ruled 
out by other observations. The ACDM model is allowed 
for p = 2 at 77% C.L. (with probability of 23%) by SN 
observations and by 0.9% C.L. (V = 0.991) by WMAP 
observations. 

In the right panel, we have shown the ratio of dark 
energy density at z = 1 and the present value. The vari- 
ation allowed by SN observations is very large, whereas 
WMAP limits the variation to within a factor 2.5 at 95% 
confidence limits. This drives the joint analysis to re- 
strict variation even further. That WMAP observations 



provide a much tighter constraint on the equation of state 
as compared to SN observations was earlier shown in |23| . 

In Fig. we show the allowed range of variation of dark 
energy as function of redshift for w = constant mod- 
els, with and without perturbations at 68%, 95% and 
99% confidence levels. The figure shows the disparity in 
allowed range by SN observations and WMAP observa- 
tions at 68% confidence level. Allowing perturbations in 
dark energy gives a similar range as compared to the case 
where dark energy perturbations are absent. In Fig. 1101 
we plot this range for varying dark energy models, top 
panel for models with p — 1 and lower panel for p = 2. 
As mentioned earlier (see also [23j), SN observations al- 
low a much wider range in change of dark energy density 
with redshift. The variation allowed by WMAP is smaller 
in all cases except constant w. The combination of the 
three constraints allows very little variation, with maxi- 
mum allowed variation in dark energy density being by 
a factor 5 up to z = 2 at 68% confidence limit. The al- 
lowed variation in dark energy density is similar in both 
the cases, indicating that the constraints on this quantity 



16 




"0 0.5 1 1.5 2^0 0.5 1 1.5 2 ~ 0.5 1 1.5 2 ~0 0.5 1 1.5 2 



FIG. 9: In this figure we plot evolution of dark energy density as a function of redshift. The top panel is for models without 
dark energy perturbations and the lower panel is with dark energy perturbations included. The left plot in both the rows shows 
the variation in energy density allowed by Supernova observations. The green/hatched region is excluded at 68% confidence 
level, red/cross-hatched at 95% confidence level and blue/solid at 99% confidence level. The white region is the allowed region 
in variation at 68% significance level. The plot which is second from left (in both the cases) displays the allowed range by 
WMAP data alone, color scheme being the same as in left figure. The next figure shows allowed range by WMAP and cluster 
abundance observations. The last column shows contribution from combined analysis of the three observations. 




~0 0.5 1 1.5 2 ^0 0.5 1 1.5 2 ^0 0.5 1 1.5 2 ^0 0.5 1 1.5 2 



FIG. 10: Same as figure 9 but for varying dark energy models. The top panel is for p = 1 and the lower panel is for p — 2. 
Starting from left, the figures show allowed dark energy density variation from SN observations, WMAP observations, WMAP 
observations combined with cluster abundance requirements and by combined analysis. The region enclosed by black solid lines 
(lower panel, second frame from left) is obtained from the derived constraints on w e ff, see text for details. 
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are parameterization independent to a large extent. 

Finally, we would like to make some comments regard- 
ing the fact that WMAP constrains the evolution of dark 
energy more effectively than SN. This arises essentially 
from the constraint on the angular diameter distance to 
the last scattering surface, or — equivalcntly the effec- 
tive equation of state parameter w e ff. (The Integrated 
Sachs- Wolfe effect and the contribution of other param- 
eters turns out to be less important.) To illustrate this 
point, we have compared the constraints on dark energy 
density for p = 2 with those implied by constraints on 
w e f f . The second figure in the lower panel in Fig. ITUI 
shows the allowed range in evolution of dark energy den- 
sity allowed by WMAP data alone for p — 2. We have 
also plotted the allowed range for dark energy density as 
a function of redshift if —1.6 < w e tf < —0.6, by thick 
black lines. This is derived by using all wq and w' that 
lead to w e ff in the range given above, and computing 
the highest and lowest dark energy density amongst this 
set of models at each redshift. We allow w and w' to 
vary in the range specified in the priors. The region al- 
lowed by the range in u» e // and that directly obtained 
from all allowed models is similar, with the latter allow- 
ing larger variation for phantom models. This reiterates 
our claim that the main constraint from WMAP data on 
dark energy parameters is on the value of W e ff at the 
last scattering surface. We believe that the larger range 
allowed at 95% confidence limit is due mainly to the ISW 
effect. 



V. CONCLUSIONS 

In this paper we presented a detailed analysis of con- 
straints on cosmological parameters from different obser- 
vations. In particular we focussed on constraints on dark 
energy equation of state, its present value and the allowed 
range of variation in it. 

It is demonstrated that the allowed range for the equa- 
tion of state parameter w is smaller if dark energy is al- 
lowed to cluster. Including perturbations mainly affects 
models with w > — 1. 

We find that WMAP observations do not distinguish 
between the ACDM model, models with a constant equa- 
tions of state parameter w and models with a variable w, 
the change in \w f° r Des t n ^ models is less than 3 even as 
the number of parameters is increased by three. WMAP 
allows only a modest variation in energy density of dark 
energy, with maximum variation being less than a fac- 
tor of three in 99% confidence limit up to z = 1. We 
infer that the main constraint from WMAP observations 
is for the derived quantity w e f / , essentially representing 
the distance to the last scattering surface. 

SN observations favour models with w < — 1 and 
O nr > 0.4. A corollory is that if we restrict to models 
with w > — 1 then the ACDM model is the most favoured 
model. Without this restriction the ACDM model is al- 
lowed only marginally by the combination of observations 



used here, this is driven mainly by SN observations. 

Allowing variation in dark energy has an impact on 
other cosmological parameters as the allowed range for 
many of these parameters becomes larger. Conversely, 
better measurements of these parameters will allow us to 
constrain models of dark energy. 

We find significant tension between different observa- 
tions. Our key conclusions in this regard may be sum- 
marised as follows: 

• SN observations favor models with large £Inr and 
w <C —1. Indeed, the best fit model is at the edge 
of our priors. 

• Enlarging priors to 0.1 < Qnr < 0.6, —0.3 > w > 
—3.0 does not lead to a better fit model for SN 
observations, indicating that our default priors are 
sufficiently wide for joint estimation of parameters. 
This is because w <C — 1 is rejected by WMAP 
observations. 

• WMAP observations favor models with w ~ — 1 
with a marginal preference for w > — 1. Including 
perturbations in dark energy removes this marginal 
preference as well. 

• For constant w models and models with variable 
w, the best fit model of each observation is ruled 
out by the other observations at a high significance 
level. As an example, the model that best fits 
the WMAP observations is completely ruled out 
by SN observations (A%| = 53). The problem is 
slightly less serious if perturbations in dark energy 
are taken into account (A%| = 12.5). 

• There is overlap of allowed regions at 95% (or bet- 
ter) by these observations, though there is little 
overlap of allowed regions at 68% confidence limit 
(see table 3). It can, of course, be argued that situ- 
ation is not alarming given that there is an overlap 
of allowed regions in parameter space at 95%. But 
we find this offset noteworthy. 

• Using larger values for erg, as indicated by some 
recent analyses ^f||6(| favors models with a slightly 
larger Qnr and slightly lower w e f / . 

• Given that the preference of individual observations 
for different types of models is not understood, and 
the fact that the best fit model of one is ruled out 
by the other, it is necessary to use a combination 
of observations for reliable constraints on models of 
dark energy. Use of either one of the observations 
is likely to mislead. 

• Our conclusions are not sensitive to priors used for 
parameters other than Qnr. Limiting priors for 
matter density to 0.1 < £Inr < 0.3 enhances the 
overlap between SN and WMAP observations and 
removes the tension between SN and WMAP ob- 
servations for constant w models. 
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TABLE III: This table lists the range of relevant parameters allowed within 68% confidence limit from SN, WMAP, WMAP + 
cluster abundance (CA) requirements and after combining all observations. 



Parameter 


AOJJ1VL 


w=const . 


w= const . 

W illi pCi I III UcLblUllb 


P-l 


p— 2 






0.256—0.3 


0.37—0.47 


0.37—0.47 


0.19—0.47 


0.18—0.47 


SN 




0.26—0.42 


0.25—0.37 


0.28—0.42 


0.22—0.39 


0.21—0.37 


WMAP 


^NR 


0.265—0.424 


0.26—0.38 


0.29—0.43 


0.24—0.42 


0.24—0.4 


WMAP+CA 




n 9K n 3 






n 9R n ^ 


0.3 0.37 


SN+WMAP+HA 






-1.9 1.39 


-1.9 1.39 


-1.89 1.31 


-1.9 1.33 


SN 






-1.0 0.67 


-1.2 0.77 


-1.19 0.56 


-1.67 0.7 


WMAP 


Wo 




-1.02 0.67 


-1.13 0.78 


-1.28 0.56 


-1.71 0.75 


WMAP+CA 






-1.41 1.02 


-1.37 1.04 


-1.3 1.05 


-1.69 1.4 


SN+WMAP+CA 










-3.2 — 3.46 


-3.3 — 3.99 


SN 










-0.94 — 0.67 


-0.74 — 4.15 


WMAP 


w'(z = 0) 








-1.0 — 0.75 
0.21 — 1.06 


-1.03 — 3.94 
2.74 — 4.58 


WMAP+CA 
SN+WMAP+CA 



TABLE IV: This table lists the range of parameters allowed within 95% confidence limit from SN, WMAP, WMAP + cluster 
abundance (CA) requirements and after combining all observations. 



Parameter 


ACDM 


w=const. 


w=const. 


P= 


1 


P= 


2 












with perturbations 














0.23 — 0.33 


0.27 - 


0.49 


0.27 — 0.495 


0.11 - 


0.49 


0.11 - 


0.49 


SN 




0.20 — 0.45 


0.16 - 


0.43 


0.20 — 0.47 


0.17 - 


0.45 


0.18 - 


0.44 


WMAP 




0.22 — 0.44 


0.21 - 


0.42 


0.22 — 0.46 


0.21 - 


0.46 


0.2 - 


0.44 


WMAP+CA 




0.23 — 0.33 


0.22 - 


0.4 


0.22 — 0.4 


0.21 - 


0.403 


0.23 - 


0.39 


SN+WMAP+CA 






-1.97 - 


-1.1 


-1.97 1.1 


-1.97- 


-1.03 


-1.97 - 


- -1.05 


SN 






-1.39 - 


-0.58 


-1.6 0.63 


-1.64 - 


-0.42 


-1.93 - 


-0.43 


WMAP 


wo 




-1.34 - 


-0.63 


-1.4 0.66 


-1.73 - 


-0.42 


-1.95 - 


-0.47 


WMAP+CA 






-1.47 - 


-0.91 


-1.48 0.96 


-1.66 - 


-0.97 


-1.89 - 


- -1.17 


SN+WMAP+CA 












-4.72 - 


- 4.6 


-4.73 - 


- 4.85 


SN 












-3.09 - 


- 1.32 


-2.5 - 


- 4.87 


WMAP 


w'(z = 0) 










-3 - 
-0.38 - 


1.35 
- 1.54 


-2.7 - 
1.15 - 


- 4.8 
4.99 


WMAP+CA 
SN+WMAP+CA 


w e ff(atLSS) 




-1.4 - 


-0.5 




-1.48 - 


-0.59 


-1.41 - 


-0.58 






0.02 — 0.027 


0.021 - 


0.028 


0.02 — 0.027 


0.02 - 


0.027 


0.02 - 


0.027 


WMAP 


Q b h 2 


0.021 — 0.027 


0.022 - 


0.028 


0.022 — 0.027 


0.021 - 


0.027 


0.022 - 


- 0.028 


WMAP+CA 




0.022 — 0.027 


0.022 - 


0.027 


0.022 — 0.027 


0.021 - 


0.027 


0.022 - 


0.028 


SN+WMAP+CA 




0.61 — 0.79 


0.61 - 


0.78 


0.6 — 0.79 


0.6 - 


0.78 


0.61 - 


0.78 


WMAP 


h 


0.61 - 0.78 


0.61 - 


0.77 


0.61 — 0.78 


0.61 - 


0.78 


0.6 - 


0.78 


WMAP+CA 




0.67 — 0.76 


0.67 - 


0.77 


0.66 — 0.76 


0.65 - 


0.75 


0.64 - 


0.72 


SN+WMAP+CA 




0.002 — 0.33 


0.011 - 


0.39 


0.007 — 0.35 


0.13 - 


- 0.4 


0.016 - 


- 0.39 


WMAP 


T 


0.001 — 0.32 


0.028 - 


-0.36 


0.006 — 0.35 


0.01 - 


0.37 


0.015 - 


- 0.387 


WMAP+CA 




0.003 — 0.3 


0.002 - 


0.29 


0.022 — 0.26 


0.006 - 


- 0.38 


0.01 - 


0.39 


SN+WMAP+CA 




0.93 — 1.08 


0.93 - 


1.1 


0.93 — 1.09 


0.93 - 


- 1.1 


0.93 - 


1.09 


WMAP 


n 


0.94 — 1.09 


0.94 - 


1.1 


0.94 — 1.09 


0.94 - 


1.09 


0.94 - 


1.09 


WMAP+CA 




0.96 — 1.09 


0.94 - 


1.08 


0.95 — 1.08 


0.94 - 


1.09 


0.93 - 


1.09 


SN+WMAP+CA 



• If we repeat the analysis of models with variable w 
with this restricted priors then we find that SN ob- 
servations strongly favor a variable w as compared 
to constant w. There is no significant tension be- 
tween SN and WMAP models in this class of mod- 
els even for the wider priors, but this is mainly due 
to a larger number of parameters. 



• Using only the Gold data set for supernovae instead 
of the Gold + Silver used here reduces the tension 
between the WMAP and supernova observations by 
a marginal amount. 

Given the points noted here regarding tension between 
fferent observations, it is important that some effort 
made to look for systematic effects in observations 



19 



as well as in analysis of observations. We have tested 
our analysis for systematic effects by varying priors and 
our findings appear to be independent of the chosen pri- 
ors, the only instance of change in results is mentioned 
above. Since the SN data set which is used by most peo- 
ple (including in this work) arises from different sources, 
one needs to be careful reg arding hidden systematics (see 
e.g., the discussion in [72|). When larger, homogeneous 
SN datasets are available in future (like for example, from 
SNLS), it is likely that the tension between the SN obser- 
vations and WMAP results disappear. If it does not, and 
the agreement continues to exist only at 3-sigma level, 



there is some cause for concern. 
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